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Abstract 

We present a theoretical study on formation rates of r\ and u 
meson- nucleus systems induced by (j,p) reactions on nuclear targets 
at ideal recoilless condition. We find that the smaller distortion effect 
in the (^y,p) reaction enables us to investigate properties of the mesons 
created deeply inside nucleus more clearly. We also consider excitation 
of scalar-isoscalar (a) mode in nucleus in order to investigate spectral 
enhancement around two-pion threshold caused by partial restora- 
tion of chiral symmetry. We conclude that valuable information of 
meson-nucleus interactions can be extracted from global structure of 
the missing mass spectra in the (7,7?) reaction. 

1 Introduction 

The study of the in-medium properties of hadrons is one of the most 
important subjects in contemporary nuclear physics and has attracted 
continuous attention. The detailed investigation of hadron- nucleus 
bound systems clarifies quantitative information on basic hadron-nucleus 
interactions. This is one of the important steps for understanding the 
QCD phase structure. So far, atomic states of pion, kaon and p have 
been investigated comprehensively and observed experimentally [1]. 
The recent interest reaches the extension to the systems with other 
heavier neutral mesons, such as the r\ and u mesons, governed purely 
by strong interaction in contrast to the atomic states. The scalar- 
isoscalar (a) excited mode in nucleus is also considered as one of the 
interesting systems. 
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Recent development in experimental aspect is the establishment 
of the (d, 3 He) spectroscopies in the formation of the deeply bound 
pionic atoms [2-5]. It opens new possibilities of the formations of 
other hadron-nucleus bound systems [6-9]. The (7,p) reaction is 
another experimental tool for formation of the hadron-nucleus system 
and originally proposed in Refs. [10-13]. The (7,p) reaction, as well as 
the (d, 3 He) reaction, is able to satisfy recoilless condition at the meson 
creation energy. The use of the (7,p) reaction has good advantage that 
the distortion effects of the projectile and the ejectile are expected to 
be smaller than in the (d, 3 He) reaction mainly due to the incident 
photon. 

It is also important mentioning that, in the theoretical point of 
view, the study of the hadron in finite density gives us information 
on the QCD phase structure [14], such as partial restoration of chiral 
symmetry, in which a reduction of the chiral condensate takes place 
in nuclear density. In this line hadronic bound systems have been 
investigated in various chiral models [6,7,15-20]. 

In this paper, motivated by the successful study of the (d, 3 He) 
reaction and the advantage of the smaller distortion effect, we discuss 
the (7,p) reaction for studying the formations of the various meson- 
nucleus systems (rj, to and scalar-isoscalar (a) mode) at the recoilless 
kinematics. Structure of the mesic nucleus is studied theoretically us- 
ing the appropriate meson-nucleus optical potentials in nuclear mat- 
ter. In order to evaluate the formation rates of the mesic nuclei, we 
introduce a Green function constructed by the Klein-Gordon equation 
with the optical potential, and calculate the cross section of the (7,p) 
reaction on a nuclear target following the method by Morimatsu and 
Yazaki [21]. In the case of the scalar-isoscalar (a) mode, we evaluate 
an effective spectral function of this mode in nucleus because of both 
the lack of the elementary cross section and the large imaginary part 
of sigma meson self-energy in vacuum, which requires certain modifi- 
cations in the Green function method. 

In the previous paper [6, 7] for studying the 77-nucleus system in 
the chiral doublet model, we found repulsive nature of the r] optical 
potential inside the nucleus which is associated with reduction of the 
mass difference of N and 7V(1535) caused by partial restoration of 
chiral symmetry. Consequently there exist certain discrepancies be- 
tween the spectra obtained by the chiral doublet model and the chiral 
unitary model. The discrepancies are expected to be distinguished by 
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experimental data [7]. However, in the (d, 3 He) reaction, due to the 
large distortion effect, it might be difficult to obtain clear information 
about the optical potential in the nuclear center. Therefore, it is in- 
teresting to consider the (7,p) reaction with much smaller distortion 
effects to explore the behavior of the i] meson in the nucleus. 

As another meson-nucleus systems, we consider the w-mesic nuclei. 
The u mesic nuclei were studied by several groups and the theoretical 
results have been reported in Refs. [8,9,12,16,17]. Recently, several ex- 
perimental results have been also reported in Refs. [22,23]. Especially, 
in Refs. [12,16,17], one finds the uo self-energy in nuclear medium based 
on an effective model and the calculated (7,p) spectra for the produc- 
tion of the (j-mesic nuclei at the recoilless kinematics. In the present 
paper, we calculate the (7,p) spectra for the formation of the w-mesic 
nuclei, and we show the incident 7 energy dependence of the whole 
spectra and the dominant subcomponents to study the experimental 
feasibilities in the lower energy photon facilities. This study is also in- 
teresting for the higher energy photon facility like SPring-8, since the 
incident photon energy has a certain distribution and we need to use 
the photon with various energies around the ideal recoilless kinematics 
to get better statistics [24]. We also study the sensitivity of the ex- 
pected spectra to the w-nucleus interaction using another theoretical 
prediction [25]. 

We also consider the (7,p) reactions in the scalar-isoscalar (a) 
channel inspired by Ref. [26-28], where an enhanced spectral func- 
tion near the 2ir threshold has been pointed out as a characteristic 
signal of the partial restoration of chiral symmetry. This is associ- 
ated with a conceivable reduction of the sigma mass and width in 
medium. Therefore, there might be a chance to create the a mesic 
nuclei with relatively small widths, if deeply bound states are formed 
in heavy nuclei [18]. It has been also found difficult, in the (d,t) or 
(d, 3 He) reactions, to observe the enhancement near the 2ir threshold 
due to the large distortion effects [18]. In this paper, we explore the 
possibility of production of the u-mesic nuclei by the photon-induced 
reactions, expecting that the 'transparency' of the (7,p) reactions en- 
ables us to observe the enhancement of the spectral function near the 
27T threshold caused by the partial restoration of chiral symmetry in 
the medium. 

This paper is organized as follows. In Sec. 2, we will discuss our 
model for the optical potential of each meson separately, and study 
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the structure of the mesic nuclei showing the binding energies and the 
widths obtained with the optical potential. In Sec. 3, we will show 
the advantages of the (7,p) reaction for the formation of the meson- 
nucleus system, and review briefly the formulation to calculate the 
formation spectrum of the mesic nucleus. In Sec. 4, we will present the 
numerical results of the missing mass spectra of the (7,p) reaction, 
and compare them with those of the (d, 3 He) reaction. Finally Sec. 5 
is devoted to summary of this paper. 



2 Structure of mesic nuclei 

In this section, we describe the theoretical formulation to study the 
structure of the meson-nucleus systems and show their numerical re- 
sults, such as the binding energies and widths, within an optical po- 
tential approach, in which all the meson-nucleus interactions are sum- 
marized in an optical potential of the meson in nuclear matter. 

The wave functions, the binding energies and their widths for the 
in-medium meson (if it is bound) are calculated by solving the Klein- 
Gordon equation with the optical potential. The Klein-Gordon equa- 
tion is written as 

_V 2 + + 2»V(u, p(r))] <f>(r) = E 2 <P(r), (1) 

where fi denotes the meson-nucleus reduced mass, which is very close 
to the meson mass for heavy nuclei, u is the real part of the relativistic 
meson energy E, and V(uj,p(r)) presents the optical potential for the 
meson-nucleus system. Here we assume to neglect the momentum de- 
pendence of the optical potential, since we consider the recoilless con- 
dition on the (7,p) reaction for the mesic nuclei formation, in which 
the meson will be created in the nucleus nearly at rest. The optical 
potential V and the eigenenergy E are expressed as complex numbers 
in general because of the absorptive effects of the nucleus. We need 
to solve the Klein-Gordon equation (1) in a self-consistent manner for 
the energy in case the optical potential V has the energy dependence. 
We follow the method of Kwon and Tabakin to solve the Klein-Gordon 
equation [29] , which was successfully applied to calculate pionic atom 
states [30]. Here, we increase the number of mesh points in the mo- 
mentum space about 10 times larger than the original work Ref. [29] 
as discussed in Ref. [6]. 
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In the following parts of the present section, we discuss the in- 
medium optical potential V(ui, p) for the rj, uo and a mesons separately. 
We show the bound state spectra, if the mesons are bound, solving the 
Klein-Gordon equation (1) with the optical potential. Here we assume 
the local density approximation with the nuclear density distribution 
p(r) written as the empirical Woods-Saxon form: 



Po 

l + exp(^) 



P(r) = ; , "r-R, * ( 2 ) 



with the radius of nucleus R = 1.18A 1 / 3 — 0.48 fm and the diffuseness 
a = 0.5 fm for the nuclear mass number A. 



2.1 77-nucleus optical potential 

First of all, we discuss the 77-nucleus system. The 77-mesic nuclei were 
studied by Haider and Liu [31] and by Chiang, Oset and Liu [32]. 
As for the formation reaction, the attempt to find the bound states 
by the (7r + ,p) reaction led to a negative result [33]. Recently, some 
experiments in photoproduction processes indicated observations of 
such bound states in 12 C target [34] and 3 He target [35], and another 
experiment is now planed to investigate the 77- nucleus interaction [9]. 
In this study we use the same theoretical models for 77-nucleus inter- 
action as described in Refs. [6, 7] in further detail. In the ry-nucleon 
system, the iV(1535) resonance (N*) plays an important role due to 
the dominant r/NN* coupling. Here we evaluate the r/-nucleus opti- 
cal potential V v (uj, p(r)) in the two different models which are based 
on distinct physical pictures of N*. One is the chiral doublet model. 
This is an extension of the linear sigma model for the nucleon and its 
chiral partner [36-38]. The other is the chiral unitary model, in which 
N* is dynamically generated in the coupled channel meson-baryon 
scattering [15,39]. 

In the first approach, the iV* is introduced as a particle with a 
large width and appears in an effective Lagrangian together with the 
nucleon field. Assuming iV*-hole excitation induced by the r\ meson 
in nucleus, we obtain the 77-nucleus optical potential at finite nuclear 
density as, 

9r/ p{ r ) 
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in the local density approximation and the heavy baryon limit [32]. 
Here \i is the r/-nucleus reduced mass and p(r) is the density distribu- 
tion of the nucleus. The r]NN* coupling is assumed to be 5-wave: 

C v nn*(x) =g v N(x)7 1 (x)N*(x) +H.c, (4) 

and the coupling constant g v is determined to be g v — 2.0 in order 
to reproduce the partial width Tn*^n ~ 75MeV at tree level. The 
S-wave nature of the i]NN* vertex simplifies the particle-hole loop 
integral in Eq. (3). m* N and m* N * are the effective masses of N and 
N* in the nuclear medium, respectively. Considering that the N* 
mass in free space lies only 50 MeV above the threshold and that 
the mass difference of N and ./V* might change in the medium, the 
77-nucleus optical potential is expected to be extremely sensitive to 
the in-medium properties of N and N*. For instance, if the mass 
difference reduces in the nuclear medium as m T] + m* N — m* N » > 0, 
then the optical potential turns to be repulsive [6]. 

Using the free space values of the N and N* masses in Eq. (3), we 
find that the optical potential gives an attractive r]N scattering length, 
a v = 0.24 + i0.38 fm, which is comparable to = 0.20-M0.26 fm and 
a v = 0.26 + i0.25 fm as reported in Ref. [39,40], respectively. On the 
other hand, the value of the r/N scattering length a v of our model is 
smaller than the recent analyses, in which larger values are reported 
such as ~ 0.9 + i0.3 fm [41,42]. In the present study, however, 
our main purpose is to study the effects of the partial restoration 
of chiral symmetry in medium. As we shall see below, in the chiral 
doublet scenario of and N* , the repulsive 77- nucleus interaction can 
be realized in the nuclear medium due to the symmetry restoration. 
This exotic behavior of the 77-nucleus interaction is independent of the 
absolute value of a v in our model. Thus, we start with the optical 
potential (3) and take into account the effects of the chiral symmetry 
restoration properly to N* and the 77-mesic nuclei. 

The chiral doublet model is used for the calculation of the in- 
medium masses of N and TV* and the in- medium iV* width in Eq. (3). 
The Lagrangian of the chiral doublet model (the mirror model) is 
given as [36, 37], 

i=l,2 

-m (iVi75iV2 - -/V275./V1) + C 

meson (5) 
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where the nucleon fields Ni and N2 have positive and negative parities, 
respectively, and the physical nucleons N and N* are expressed as a 
linear combination of Aq and N2. The parameters are determined 
so that the model reproduces the free space values of the N and iV* 
masses and the partial decay width of N* to irN as Tn*->wN = 75 
MeV with (<r) = 93 MeV [37]. 

In the chiral doublet model together with the assumption of partial 
restoration of chiral symmetry, a reduction of the mass difference of N 
and N* in the medium is found to be in the mean field approximation 
as [36,37,43,44], 

m* N (p) - m* N *{p) = $(p)(m N - m N *), (6) 

where win and tun* are the N and N* masses in free space, respec- 
tively, and 

Hp) = 1 _ C ^. (7) 
Po 

Here we take the linear density approximation of the in-medium mod- 
ification of the chiral condensate, and the parameter C represents the 
strength of the chiral restoration at the nuclear saturation density po- 
The empirical value of C lies from 0.1 to 0.3 [27]. Here we perform 
our calculations with C = 0.0 and 0.2 in order to investigate the ef- 
fect of the partial restoration of chiral symmetry. A self-consistent 
calculation within the chiral doublet model gives moderately linear 
dependence of the chiral condensate to the density with C = 0.22 in 
Ref. [44], in which the relativistic Hartree approximation is used to 
calculate the chiral condensate in nuclear matter. 

In Fig. 1, we show the rj- nucleus optical potentials obtained by 
the chiral doublet model with C = 0.0 and C = 0.2 in the case of 
the r/- 11 B system assuming the local density approximation and the 
Woods-Saxon nuclear distribution defined in Eq. (2). The potential 
with C = 0.0 corresponds to the so-called tp approximation because 
there are no medium modifications for N and iV*. As shown in the 
figure, the potential with C = 0.2 has quite different shape from the 
case with C = 0.0. The r/-nucleus optical potential with C = 0.2 
turns to be repulsive above a critical density p c at which the sign of 
lo + m* N — m* N * becomes positive. Consequently, the 77-nucleus optical 
potential has a curious shape of a repulsive core inside nucleus and 
an attractive pocket in nuclear surface as indicated in Fig.l. The 
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qualitative feature of the optical potential discussed here, such as the 
appearance of the repulsive core in the case of C = 0.2, is independent 
of type of nuclei due to saturation of nuclear density. It is extremely 
interesting to confirm the existence (or non-existence) of this curious 
shaped potential experimentally. 

We should mention here that there are two types of the chiral dou- 
blet model, namely, the mirror assignment and the naive assignment, 
due to two possible type of assignment for the axial charge to the iV 
and N* [36,37]. The potential in Fig. 1 is obtained by the mirror as- 
signment case. In the mirror assignment, N* is regarded as the chiral 
partner of N and forms a chiral multiplet together with N. While in 
the naive assignment, N and N* are treated as independent baryons 
in the sense of the chiral group. We calculate the reaction spectra 
using both assignments, and show the results in next section. In both 
cases, the mass difference of N and N* is given in the same form 
as shown in Eq. (6). Hence, the energy and the density dependence 
of the potential of the naive assignment resemble those of the mirror 
assignment. The detail discussions are given in Ref. [7]. 

Let us move on the second approach of the chiral unitary model [15, 
39]. In this approach, the N* resonance is expressed as a dynami- 
cally generated object in the meson baryon scattering, and one solves 
a coupled channel Bethe-Salpeter equation to obtain the ?7-nucleon 
scattering amplitude. The optical potential in medium is obtained by 
closing the nucleon external lines in the r/N scattering amplitude and 
considering the in-medium effect on the scattering amplitude, such as 
Pauli blocking. Since the N* in the chiral unitary approach is found 
to have a large component of KT, and the S hyperon is free from the 
Pauli blocking in nuclear medium, very little mass shift of N* is ex- 
pected in the medium [45] , while the chiral doublet model predicts the 
significant mass reduction as discussed above. The optical potential 
obtained by the chiral unitary approach in Ref. [15] is also shown in 
Fig. 1, and the potential resembles that of the chiral doublet model 
with C = 0.0 due to the small medium effect. 

In Ref. [6], in order to obtain the eigenenergies, we have solved the 
Klein-Gordon equation Eq. (1) numerically with the 77-nucleus optical 
potential given in Eq. (3) for the mirror assignment. We show in 
Table 1 the calculated binding energies and level widths for 0s and 
Op states in the 11 B and 39 K cases. We also show in Table 1 the 
binding energies and widths by the chiral unitary approach reported 
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Figure 1: The 77-nucleus optical potential for the r/- n B system as functions 
of the radius coordinate r reported in Ref. [7]. The left and right figures 
show the real part and the imaginary part of the ^-nucleus optical potential, 
respectively. In both figures, the solid lines show the potentials of the chiral 
doublet model with C = 0.2 (thick line) and C = 0.0 (thin line), and dashed 
lines show that of the chiral unitary model, which is picked from the results 
shown in Ref. [15] 



in Ref. [15] for comparison. The level structure of the bound states is 
qualitatively very similar with that of the chiral doublet model with 
C = 0.0 as we can expect from the potential shape. For the C = 0.2 
case in the chiral doublet model, we can not find any bound states 
because of the repulsive nature of the optical potential. 

2.2 cj-nucleus optical potential 

Let us move to discussion on the optical potential of in-medium ui 
meson. In our study of the ui mesic nuclei, we use a simple model 
for the optical potential, which was used in the previous study [9] on 
the w-mesic nuclei formation using the (d, 3 He) reaction. An empirical 
potential is given in an energy-independent Tp approximation as 

V^(r) = -(100 + 70z)^. (8) 

Po 

The real and imaginary parts of the optical potential were estimated 
based on 15 % mass reduction of the u mass at the normal density 
and the ui life time (r = 1.5 fm/c) of nuclear absorption in the nuclear 
medium, respectively, as reported in Ref. [16]. 
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Chiral doublet model with C=0.0 [MeV] Chiral doublet model with C=0.2 

"B 39 K = ^ 1 B 39 K 

Os (13.7, 41.5) (30.3, 42.5) no bound state 

Op - (14.6, 50.7) 



Chiral unitary approach [15] [MeV] 



12 C 


4u Ca 


0s (9.71, 35.0) 


(17.9, 34.4) 


Op 


(7.0, 38.6) 



Table 1: Binding energies and widths calculated by the Chiral Doublet model 
(C=0.0 and C=0.2) with the mirror assignment [6] and by the chiral unitary 
approach reported in Ref. [15]. 



u- n B bound states 


T/(a) 
'id 


y(b) 
'Ui 


0s 


(49.0, 116) 


(103, 54.4) 


Is 




(18.8, 27.6) 


Op 


(14.3, 86.5) 


(58.1, 42.7) 


Od 




(17.0, 30.5) 



Table 2: Binding energies and widths of the u bound states in n B with the 
optical potentials (8) and (9) in unit of MeV. 

The binding energies and widths of the omega meson in n B are 
calculated by solving the Klein-Gordon equation (1) with the above 
potential. In this potential, we find two bound states in the n B nu- 
cleus. The obtained binding energies and widths are shown in Table 2. 
Since the width of these states are larger than the level spacing, it 
would be difficult to observe distinct peak structure in inclusive spec- 
tra. In next section, we will show the calculated (7,p) spectra for 
several incident energies. 

We also consider two more optical potentials in order to investigate 
the sensitivity of the depth and the absorption of the potential to the 
(t>p) spectra. The potentials considered here are 

vS°\r) = -(156 + 29i)^ , (9) 

Po 
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yj c )(r) = _(_42.8 + 19.5i)^ (10) 

Po 

which are obtained by the linear density approximation with the scat- 
tering lengths a = 1.6 + 0.3i fm [17] and a = -0.44 + 0.2i fm [25], 
respectively. The former scattering length is obtained by the loN 
scattering amplitude calculated at one loop approximation with an ef- 
fective Lagrangian based on the SU(3) chiral symmetry incorporating 
vector mesons. On the other hand, the later is obtained in a relativis- 
tic and unitary approach for meson-baryon amplitudes. In this case 
the scattering length is to be repulsive due to a subthreshold effect of 
the wiV(1520). 

The energies and widths of the bound states calculated with the 
potential vS°^ are also shown in Table 2. In this case, more bound 
states are found with relatively narrower widths than the case of the 
previous potential vS^ since the optical potential vS°^ has deeper real 
potential and less nuclear absorption. The magnitudes of the level 
widths are comparable to those of level spacing and, hence, we expect 
to observe peak structure in reaction spectra. Later we shall com- 
pare the (7,p) spectra calculated with these potentials at the incident 
energy E 1 = 2.7 GeV where the recoilless condition is satisfied. 



2.3 cr-nucleus optical potential 

Finally, we consider the a meson-nucleus systems. In this paper we 
take a particle picture of the sigma meson, in which the sigma meson 
is described as a particle with a huge decay width to two pions. The 
calculation of the self-energy of the sigma meson is performed based 
on the SU(2) linear sigma model [27], and the partial restoration of 
chiral symmetry in nuclear medium is parameterized as the sigma 
condensation in the medium within the model. 

The optical potential for the cr-nucleus is defined as the density- 
dependent part of the in-medium a self-energy divided by two times 
the bare sigma mass m a as; 

BeV a (u,p) = "^-Acr 2 (1 - <D(p)) . (11) 

Here we use the a self-energy evaluated in Ref. [27] based on the 
SU(2) linear sigma model at the one-loop level, which corresponds to 
the mean field approximation. The density dependent function <&(p) 
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is defined to express the chiral condensate in nuclear matter (a) p as 
(a) p = ao&(p). The function $(p) is given in Eq.(7) with the C 
paramater, which represents the strength of the partial restoration of 
the chiral symmetry in the nuclear medium. It is worth noting that 
at the one-loop approximation the imaginary part of the self-energy is 
independent of the density, so that the nuclear absorption of the sigma 
meson appears from the next order, such that multi-p/i excitations and 
pion - ph excitations. 

The density-independent part of the self-energy E CT is given in the 
one-loop calculation in the linear sigma model as 



ReY, a (u>) = - 



A 



32ir 2 



ml(l - In ^f) + ml (1-ln^) 



K 



ml 



+ ^ (Q w + 2-ln-f0 



ml 



+ iQcr + 2 — In — 



K 



Im£ CT (u,)= - A- 



4 



\p n e(uj - 2m n ) + P a e{uj - 2m a ) 



(12) 
,(13) 



where the function Q v is given as 



(for 2m v < lo 



-2P V arctan — (for uj < 2m if ) 
Pip 



(14) 



with P w 



1/2 



Here, ip denotes either tt or a, k is a 
renormalization point in the minimal subtraction scheme and A is the 
coupling constant of the four-point meson vertex in the linear sigma 
model. The reduced mass fi shown in the Klein-Gordon equation (1) 
is corresponding to fi 2 = m\ + ReSo- in heavy nucleus limit for the 
sigma meson case. The bare sigma mass m a is fitted as 469 MeV 
so that the sigma meson mass in free space mh ee ^ is reproduced as 
550 MeV after the one-loop calculation. All parameters appeared in 
Eqs. (12)-(14) are obtained in the first reference of Ref. [26] and listed 
in Table. 1 of Ref. [18]. 

The self-energy of the sigma meson obtained in the mean field pro- 
vides a strong attractive potential for the sigma meson in medium. We 
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show in Fig. 2 the attractive potential ReV a (r), given in Eq.(ll), for 
C = 0.2, 0.3 and 0.4 in 208 Pb. We assume that the density distribu- 
tion is of the Woods-Saxon type as defined in Eq. (2) with R = 6.5 fm 
and a = 0.5 fm. It is seen that the depth of the attractive potential 
is strongly depend on the C parameter. The strong attraction of the 
optical potential in the center of nucleus appears also in lighter nuclei. 

As a result of the attraction, we expect that a strong reduction of 
the a mass approaching down to the two pion threshold especially for 
heavy nuclei. Consequently there is a chance to form sigma bound 
states with narrower widths. The reduction of the sigma mass in 
nuclear medium is expected also in the irir resonance picture of the 
sigma meson with the non-linear realization of chiral symmetry, if one 
takes account of appropriate medium effects through the pion wave 
function renormalization [28]. 

The energies and widths of the sigma bound states are calculated 
by solving the Klein-Gordon equation for the a meson with the optical 
potential V a and the self-energy S CT with the one-loop correction. In 
Fig. 3, we show the bound state spectra of the a meson in 208 Pb with 
C = 0.2, 0.3 and 0.4 cases, together with widths of low lying states [18]. 
It is found that, if the partial restoration of chiral symmetry in nuclear 
medium occurs with sufficient strength, deeply bound a states in heavy 
nuclei formed with significantly smaller decay widths than that in 
vacuum. This is responsible for suppression of the available itit phase 
space associated with the reductions of the real energies of the bound 
states. Furthermore, for C = 0.3 and 0.4 cases, there are several 
bound states with no decay widths below the 2n threshold since the 
a — > 7T7T decay channel, which is the only possible decay mode at low 
energies in the present model, is kinematically forbidden. 

In order to see a global structure of the a bound states in the 
nucleus and the impact of the narrow bound states on the production 
reaction, we show the spectral functions of the bound states, which 
are used to calculate the a production spectra in the (7,p) reaction 
later. The spectral function of the (n, £) bound state is given by 

/ ^_ 1 ImS g (tj) 

Pn£iUj) ~ -^-^ + (Im£ CTH ) 2 ' (15) 

where u n £ denotes the eigenenergy of the a bound state in the nucleus. 
The spectral functions Eq. (15) encounter a divergence at u = uj h £ for 
eigenstates in case of ImS CT = 0. This is not the case if we take 
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Figure 2: Real potential for the a meson inside the 208 Pb nucleus as a function 
of the radial coordinate r obtained in Ref. [18]. The density is assumed to 
be the Woods-Saxon form with R = 6.5 fm and a = 0.5 fm in Eq. (2). Each 
line indicates C = 0.2 (dotted line), 0.3 (dot-dashed line) and 0.4 (solid line) 
case, respectively. 

into account nuclear absorptions of the sigma meson which we do not 
consider here. In order to regularize the divergence, we add an extra 
5 MeV width to the ImS ff in all the bound states in the following 
calculations. 

Figure 4 shows a series of p n i for the t = bound states in 208 Pb 
for the C = 0.3 case [18]. As seen in the figure, the spectral functions 
for the deeply bound states have a prominent peak around u ~ 2171^. 
This is a direct consequence that the deep bound states are formed 
as a result of the strong attractive potential. The detail discussions 
are given in Ref. [18]. The appearance of the narrow states at the low 
energies gives also strong enhancement of the sigma production rate 
in nucleus around the two pion threshold, which is expected to be seen 
in the (7,p) reaction. This point is discussed in the following section. 

It is also possible to form the deeply bound states with the narrow 
width in lighter nuclei than Pb, since the potential depth is nearly 
independent of the mass number of nucleus for a certain value of the 
C parameter. The spatial range of the potential, however, depends on 
nuclear size. Therefore, the zero width bound states are expected to 
be formed in nuclei which are larger than a certain spatial dimension. 
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Figure 3: Calculated eigenenergies of the bound a in 208 Pb [18] with (a) 
C = 0.2, (b) C = 0.3 and (c) C = 0.4. The £ is the orbital angular momentum 
of the a states. Widths are shown by numbers for low lying states in unit of 
MeV. The 2m n threshold is shown by the dotted line. The states blow the 
threshold do not have widths due to the mr decay in the present calculation. 

3 Formation of meson-nucleus system 
in (7,p) reaction 

In this section we explain the advantages of the (7,p) reaction for 
investigation of the meson-nucleus system and briefly review the for- 
mulation to obtain the formation rate of the mesons inside nucleus. 



3.1 (t ? p) reaction 

In the (7,p) reaction, the incident photon reacts with a proton inside 
the target nucleus, producing a meson there most likely at rest with 
recoilless condition. The meson may be trapped into a certain bound 
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Figure 4: Spectral functions p n i(uj) of the a bound states with i = and 
n = 1,2, ... ,5 for C = 0.3 case obtained in Ref. [18]. An extra 5 MeV width 
is added to the imaginary part of the self-energy. 



state in the nucleus, if there exist bound states. The proton is ejected 
from the nucleus and is to be observed at a forward angle. In the 
(7,p) reaction spectroscopies, one observes only the emitted proton 
in the final state for simplicity of the experiment and obtains the 
double differential cross section da /dfl /dE as a function of the emitted 
proton energy. Investigating the differential cross section, we extract 
the properties of the interaction between each meson and the nucleus. 
Observing the more particles in the final state, we will obtain the more 
information of the meson-nucleus interactions. 

The incident energy for each meson production is tuned so as to 
satisfy the recoilless condition achieved by the zero momentum trans- 
fer, in which the meson may be produced in the nucleus at rest. 
In the recoilless condition, due to the matching condition between 
momentum- and angular momentum-transfer, we have a selection rule 
for the total angular momentum of the meson and the proton hole 
states. Namely the contributions of the different angular momenta 
are strongly suppressed. 

In Fig. 5, we show the momentum transfers as functions of the 
incident particle energies for the (d, 3 He) and (7,p) reactions. In the 
plots we use the meson masses in free space. In the following discussion 
for the (7,7?) reaction we choose the incident energy = 950 MeV for 
the eta meson production and J3~ = 2.75 GeV for the omega meson 
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Figure 5: Momentum transfer q is shown as a function of (a) the incident 
deuteron kinetic energy T d in the (d, 3 He) reaction, and (b) the incident 7 
energy E 1 in the (7,p) reaction. The solid line shows the momentum transfer 
of the pion production case, and other lines show those of two pions, the r\ 
meson and the ui meson production cases as indicated in the figure. 

production to satisfy the recoilless condition for the corresponding 
meson-nucleus systems. For the sigma meson case, the incident energy 
is chosen to be = 400 MeV, which is a recoilless condition for 
the two pion production in the nucleus, since we focus our interest 
on the spectral function enhancement near the two pion threshold 
in the sigma channel associated with the partial restoration of chiral 
symmetry in nuclear medium. 

The major advantage of use of the (7,p) reaction for the mesic 
nucleus formation is the smaller distortion effect on the incident pho- 
ton than hadronic beams. The distortion effects of the projectile and 
ejectile are known to reduce significantly the reaction rate inside the 
nucleus. As we shall discuss later in Fig. 6, half of the flux of the (7,p) 
reaction reaches the center of the target nucleus, while most of the 
flux in the (d, 3 He) reaction is restricted on the surface of the nucleus. 
Therefore, due to the transparency of photon, we can probe the me- 
son wavefunction in the center of the nucleus with the (7,p) reaction. 
This characteristic feature is desirable especially to the study of the 
eta-nucleus system, in which we have the discrepancy in the model 
calculations of the 77 optical potential at the center of the nucleus 
depending on the physical picture of N* , as discussed in Sec. 2.1 

The other advantage of the (7,p) reaction is that, due to the trans- 
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parency of photon, we expect more production rate of the low lying 
states of the mesic nucleus. This is favorable to the sigma-nucleus sys- 
tem, since we expect the deeply bound states of the sigma meson with 
a significantly narrow width as a result of strong reduction of the chiral 
condensate in the nucleus. In the case of mesic nuclei formation, the 
meson wave function in nucleus resembles the proton hole wave func- 
tion in the sense of function orthogonality, because both of the wave 
functions are, roughly to say, solutions of the Schrodinger equations 
with the same square-well type potentials. (Recall that the potentials 
for the mesons is obtained as the Woods-Saxon type shape under the 
local density assumption with the nuclear density distribution (2) ex- 
cept the case of the eta meson with the chiral doublet model C = 0.2.) 
Due to the orthogonality of the wave functions, we have the selection 
rule that the combinations of the different principle quantum numbers 
of the meson and proton hole wave functions are strongly suppressed 
in the recoilless condition [9,18]. This is reminiscent of the recoilless 
formation of substitutional states in hypernuclear reactions. There- 
fore, in order to produce the ground state of the mesic nucleus, in 
which the meson sits on the s state, we need to pick up the deepest 
bound proton from the target nuclei using a better transparent inci- 
dent beam. In contrast, in the case of the atomic state formation, 
the orthogonal condition of the radial wave functions of the meson 
and nucleon is not satisfied any more, since the meson is trapped in a 
Coulomb potential and the meson distribution is roughly expressed by 
the Coulomb wave function. Then the recoilless condition enhances 
the population of any states which have the total angular momentum 
J ~ 0, for example, in the formation of the deeply bound pionic atoms 
in Sn isotopes [5,46], the contribution of (3s)" 1 ® ls w combination can 
be largely populated in the (d, 3 He) reaction. 

Let us evaluate the distortion effects of the (d, 3 He) and (7,p) reac- 
tions to compare the 'transparency' of both reactions. We describe the 
distorted waves of the incoming projectile and of the outgoing ejectile 
as Xi an d X/> respectively. Using the Eikonal approximation, we write 

X/( r ki(r) = exp(iq • r)F(r), (16) 

with the momentum transfer q. We choose the beam direction as to 
be z-axis. The momentum transfer q is along z-axis in the forward 
reactions which we consider here. The distortion factor F(r) is defined 
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as, 

r 1 f z 1 f°° 

F(r)=exp --(JiN dz'p A (z',b) - -a fN dz'p A -i(z', b) . 

Z J— oo Z Jz 

(17) 

where a in (c/tv) is the total cross section of the nucleoli and incident 
(emitted) particle and p A (z, b) is the density distribution function for 
the nucleus with the mass number A in the coordinates of z = r cos 9 
and b = rsin#. The distortion factor, which will be used in the 
calculation of the (7,p) reaction cross section later, gives an estimation 
of the reduction of the beam flux caused by the nuclear absorption 
effects. In order to take a view of the distortion effect, we calculate 
an averaged distortion factor, which is defined as 



F(b) = exp 



1 r°° 

-a p(z',b)dz' 

Z J— oo 



(18) 



2 

where a is an average distortion cross section of the initial and final 
channels, and p(z, b) denotes an average nuclear density of the target 
and daughter nuclei at an impact parameter b and a beam direction 
coordinate z. 

In Fig. 6, we show the averaged distortion factors F(b) in the (7,p) 
reaction and the (d, 3 He) reaction of the 12 C target case for r) produc- 
tion, which will be discussed in Sec. 4.1. Here we have used a = 104 
mb for the (d, 3 He) reaction and a = 16 mb for the (7,p) reaction. 
These cross sections are evaluated from experimental nucleon-nucleon 
scattering data at appropriate energies [47]. The figure shows that 
F(b) in the (7,p) reaction has a finite value (> 0.5) even in the nu- 
clear center. This means that the amplitude of the meson production 
at center of nucleus is suppressed only by half due to the nuclear ab- 
sorption. On the other hand, the averaged distortion factor of the 
(d, 3 He) reaction is almost zero inside the nucleus. Therefore the pho- 
ton can reach the center of the nucleus and create the meson there, 
while the deuteron mostly creates meson only on the surface of the 
nucleus. Thus, as we expected, the (7,p) spectra are shown to be more 
sensitive to the optical potential of the meson at the interior of the 
nucleus. 



3.2 Formulation 

We use the Green function method [21] to calculate the formation 
cross sections of the ry-nucleus and the cj-nucleus systems in the (7,p) 
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reaction. The details of the application of the Green function method 
for the (d, 3 He) reaction are found in Refs. [6,9, 17]. 

The present method starts with a separation of the reaction cross 
section into the nuclear response function S(E) and the elementary 
cross section of the p(7,p)?y or uj with the impulse approximation: 

da \ lab 

I x S(E) , (19) 




A( 7 ,p)(A-l)® ¥) / p(j,p)<p 

where 93 denotes the r\ or u> meson. The calculation of the nuclear 
response function with the complex potential is formulated by Mori- 
matsu and Yazaki [21] in a generic form as 

S(E) = --lmJ2T}G(E)T f (20) 
/ 

where the summation is taken over all possible final states. The am- 
plitude Tf denotes the transition of the incident particle to the proton 
hole and the outgoing ejectile, involving the proton hole wavefunction 
ipj p and the distorted waves, \i an d Xji °f the projectile and ejectile, 
taking the appropriate spin sum: 

7)(r) = X}(r)e i/2 , ms ®lMr)] JM *« (21) 

with the meson angular wavefunction Yi (f) and the spin wavefunc- 
tion £i/2,m s °f the ejectile. The distorted waves are written with 
the distortion factor as Eq.(16) in the Eikonal approximation. The 
Green function G(E) contains the meson-nucleus optical potential in 
the Hamiltonian as 

G(E;r,r') = (p-i|^( r )__±— ^Ol^ 1 } (22) 

where is the meson creation operator and is the proton hole 

state. Obtaining the Green function with the optical potential is essen- 
tially same as solving the associated Klein-Gordon equation. We can 
calculate the nuclear response function S{E) from Tj(r)G(E; r, r')T / (r') 
by performing appropriate numerical integrations for variables r and 
r'. 

For the cr-nucleus systems, we do not use the Green function 
method to evaluate the formation rate, because (a) we do not know 
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the elementary cross section of the 7 + p — > a + p reaction yet and 
(b) the sigma optical potential has a large imaginary part even out- 
side of the nucleus due to the large decay width, of which boundary 
condition the Green function does not work in the present formalism. 
Alternatively we evaluate the spectral function of the a bound states 
defined as 

p a (E) = ]T(2/ + l) Pni (E) , (23) 

n,e 

where p n e(E) denotes the contribution from each bound state in the 
nucleus with the (n,£) quantum number as given in Eq. (15). In 
order to take the reaction mechanism into consideration, we introduce 
an effective number N e g, which represents relative production wight 
of each bound state [2, 18]. Finally we obtained the total spectral 
function as 

Ptot(^) =J2N eSPne (u). (24) 

nl 

with the effective number N e g given by 

N cS = E I / d 3 r X }(r)e i/2 , ms [#» ® V^M] JM Xi(*)\ 2 (25) 

,JMm s J 

where tpj p denotes the wave function of the picked-up proton and, <^ CT 
and £i/2,m s are the a wave function and the spin wave function of the 
ejectile, respectively. The a wave function is obtained by solving the 
Klein-Gordon equation with the optical potential and the self-energy 
shown in Eqs. (11), (12) and (13). 

This effective number approach is known to be a good approxi- 
mation of the Green function method to calculate the reaction cross 
section in Eq.(19), if the bound states lie well-separately with narrow 
widths [21], and applied to evaluate the formation cross sections of 
deeply bound pionic atoms [2]. In the present case of the sigma me- 
son in nucleus, we expect the narrow bound states around the two-pion 
threshold for the C = 0.3 and C = 0.4 cases, and focus our interest on 
seeing the impact of the narrow bound states on the (7,p) reaction. 
Therefore the approach with the spectral function is expected to be 
reasonable for our purpose. 
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4 Result 



In this section, we evaluate the formation rates of the mesic nuclei by 
the (7,p) reactions and show the spectra calculated with the optical 
potentials discussed in Section 2 for the rj-, u-, and a-nucleus systems. 

4.1 77-nucleus system 

We discuss in this section the formation rate of the 77-mesic nucleus 
in the (j,p) reaction. The initial photon energy is chosen as the 
recoilless condition energy = 950 MeV as shown in Fig. 5(b). The 
emitted proton spectra around the rj production energies are evaluated 
in the Green function method (19) with the 77 Green function (22) 
calculated with the optical potential (3). We estimate the elementary 
cross section da / dVL of the 7 + p — > p + rj reaction at this energy in 
the laboratory frame as to be 3.4/ib/sr using the experimental data 
by the tagged photon beam at University of Tokyo [48] . 

The obtained spectra are shown as functions of the excitation en- 
ergy E ex measured from the eta meson production threshold energy 
Eq. The excitation energy corresponds to the missing energy in the 
inclusive reaction 7 + A — > p + X . In the plot of the separated contri- 
butions of the 77-hole configurations, it should be taken into account 
the appropriate hole energy depending on the orbit j p in the daugh- 
ter nucleus. The hole energy S p {j p ) — Sp(ground) is measured from 
the ground state of the daughter nucleus, which are listed in Table 3. 
The binding energy B v of the 77 meson is deduced from the excitation 
energy and the hole energy as 

E cx = 771^-3^ + [S p (jp) - S p (ground)]. (26) 

We evaluate the differential cross section of the (7,p) reaction on 
the 12 C target. The r\ meson is created in the 11 B nucleus because 
the proton is knocked out in the reaction. In Fig. 7, we show the 
results calculated by the chiral doublet model with both mirror and 
naive assignments together with those of the (d, 3 He) reaction reported 
in [6,7] for C = 0.2 case in Eq.(7). 

First of all we would like to mention that the peak structure seen in 
the spectra has nothing to do with the existence of the rj bound states. 
(Actually in the case of C = 0.2 there are no bound states.) These 
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Table 3: Proton separation energies of the single particle orbits measured 
from the ground state of the daughter nucleus are shown in unit of MeV. 
Data taken from Ref. [49] for 12 C and Ref. [50] for 40 Ca. Widths of the 
proton hole states deduced from the data are also shown for 40 Ca. 
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Figure 6: Distortion factor F(b) for the r) meson production reactions in 12 C 
target as a function of the impact parameter b. The solid line indicates the 
F(b) for the (7,p) reaction and the dashed line for the (d, 3 He) reaction. The 
density profile of 12 C is also shown. 

peaks are just signals of opening the phase space of the rj creation in 
the nucleus [7]. 

Secondly we see in the figure that each spectrum is dominated by 
two contributions, (Os^)" 1 (8> s v and (0p3/2) p ~ 1 ®Pt] This is a conse- 
quence of the matching condition in the recoilless kinematics that final 
states with total spin J ^ are largely suppressed. The rj production 
threshold for the {0p3/2) p ~ 1 proton-hole state is indicated as vertical 
thin line at E ex — E$ = 0. The threshold for the (0si/2)p 1 hole state, 
which is the excited state of the daughter nucleus, is at E ex — Eq = 18 
MeV because of the different S p (j p ) in Eq. (26). 

As we can see in the figure, the difference between the naive and 
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mirror assignments is enhanced in the (7,p) reaction, as compared 
with the (d, 3 He) case and, hence, the (7,p) reaction is more sensitive 
to the details of the fy-nucleus interaction as we expected. However, it 
might be difficult to distinguish these two cases by experimental data, 
since the difference is seen only in the magnitude of the spectra above 
the rj production threshold. Hereafter, we show only the results with 
the mirror assignment. 

In Fig. 8, we show again the 12 C target cases for the (7,p) and 
(d, 3 He) reactions of the r\ mesic nucleus, comparing the three differ- 
ent 7/-nucleus optical potentials. In Fig. 8(a), the spectra with the so- 
called tp optical potential are shown, which are calculated by putting 
C = 0.0 in the chiral doublet model. We show also the spectra ob- 
tained by the chiral doublet model with C = 0.2 in Fig. 8(b). Compar- 
ing (a) and (b) in the left-side figure, we find that the repulsive nature 
of the potential in the chiral doublet model with C = 0.2 makes the 
bump structure of the (d, 3 He) spectrum broaden out to the higher en- 
ergy region, as discussed in Ref. [6]. On the right-hand side in Fig. 8, 
the same 'broadening effect' is seen also in the (7,p) reaction. We 
would like to stress here that this tendency to widen the spectra is 
seen more clearly in the (7,p) spectra as a result of the transparency 
of the incident photon. In Fig. 8(c), we show the spectra calculated 
by the chiral unitary approach. We can see that the spectra of the 
chiral unitary approach are shifted to the lower energy region like in 
Figs. 8( cl^ clS cl result of the attractive potential. It is very interesting 
to see which of spectra, (b) or (c), is realized in experiment. It is 
also interesting to compare the spectra of the (7,p) and (d, 3 He) reac- 
tions in experimental observation, since we can expect to obtain the 
detail information of the r/-nucleus interaction by the comparison of 
the both data which have different sensitivities to the nuclear surface 
and center. 

As for the contributions from the bound r\ states in the (7,p) re- 
action, in Fig. 8(a) and (c), we find a certain enhancement in the 
(0si/2)p 1 ® s v configuration as a bump structure below the 0s state 
eta production threshold, which is E cx — Eq = 18 MeV. This is the 
indication of the rj meson bound in the 0s state, which is listed in 
Table 1. Although the bound state is more clearly seen in the (7,p) 
reaction, in which we expect more reaction rates to pick up a proton 
in deeper states in the target nucleus, as discussed before, it is still 
hard to distinguish the bound state as a peak structure even in the 
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(7,p) reaction due to the large width of the bound state and the larger 
contribution of (0p3/2)p 1 state. 

Next we consider the heavier target 40 Ca for the formation of the r\ 
mesic nucleus. In this case, we have larger possibilities to have bound 
states and to observe larger medium effects. In Fig. 9, we show the 
calculated spectra of 40 Ca(7,p) 39 K<g)77 reaction using the chiral doublet 
model with C = 0.0 and C = 0.2 together with those of the (d, 3 He) 
reaction for comparison. As in the case of the 12 C target, we find larger 
difference between the spectra with C = 0.0 and C = 0.2 in the (7,p) 
reaction than those in the (d, 3 He) reaction. In the (d, 3 He) spectra, 
the total spectra are shifted around 5 MeV to higher excitation energy 
region with almost same shape by changing the C parameter value 
from C = 0.0 to C = 0.2. On the other hand, the (7,p) spectra 
are largely shifted to higher energy region with significantly deformed 
shape. This is the good advantage of the (7,p) reaction to deduce the 
information of the rj optical potential. We can also see the contribution 
from the bound r) states for the C = 0.0 case more clearly in the (7,p) 
spectra. The bound states were obtained in the C = 0.0 case as energy 
30.3 MeV with width 42.5 in s-wave and energy 14.6 with width 50.7 
in p-wave [6]. However, as in the case of the (d, 3 He) reaction, it seems 
difficult to deduce the contribution from a certain subcomponent from 
the total spectra even in the (7,p) reaction because of the large widths 
of the r? states. 

In short summary for the 77 mesic nucleus formation, we find that 
there exist certain discrepancies between the spectra obtained with 
the different chiral models and we also find that the (7,p) spectra are 
more sensitive to these discrepancies than the (d, 3 He) spectra. 

4.2 cj-nucleus system 

In the uj production in the (7,p) reactions, the recoilless condition is 
satisfied by the incident photon energy E 1 = 2.75 GeV. The elemen- 
tary cross section of the 7 + p — ► p + uj process is evaluated to be 
0.3/jb/sr [12]. 

First of all we show the result of the ui mesic nucleus formation 
spectra with the optical potential given in Eq. (8). Shown in 
Fig. 10 are the calculated spectra of the (7,p) reaction on the 12 C 
target with various incident energies for = 1.5, 2.0, 2.7 and 3.0 
GeV. So far the formation rate was calculated for the ideal kinemat- 
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Figure 7: Calculated ^-production spectra of the (a) (d, 3 He) reaction at 
T d = 3.5 GeV and (b) (7,p) reaction at £ 7 = 950 MeV for the 12 C target 
are shown as functions of the excitation energy E ex defined in the text. E is 
the i] production threshold energy. The ^-nucleus interaction are evaluated 
using the chiral doublet model with the mirror assignment (thick lines) and 
naive assignment (thin lines) with C = 0.2. The total spectra are shown by 
the solid lines and the dominant contributions from the (0si/2)p 1 <S> s v and 
(0^3/2)p 1 ®Pr) configurations are shown by the dashed lines and dash-dotted 
lines, respectively, where the proton-hole states are indicated as (nij)' 1 and 
r\ states as i^. 



ics to satisfy the recoilless condition. Here we study also the photon 
energy dependence of the expected spectra and its subcomponent con- 
tributions so as to see the experimental feasibility in the lower energy 
photon facilities. The energy dependence of the expected spectra is 
also useful for the planned experiment at SPring-8 because it is re- 
quired to use various photon energy around the ideal kinematics to 
have the data with better statistics [24]. As shown in Fig. 10, the spec- 
tra are dominated by the contributions from the two configurations 
(0si/2)p 1 ® $u) and ifipz/2)p l ® Pu for the incident energies E 1 = 2.7 
GeV and 3.0 GeV, which roughly satisfy the recoilless condition. Even 
in this ideal case for clear observations, the bound state structure is 
hardly seen in the spectra due to the large widths of the bound states. 
At £L = 2.0 GeV, the total spectrum is still dominated by the two con- 
figurations but other subcomponents have certain contributions above 
the u production threshold. For E~ = 1.5 GeV, we find that many 
configurations contribute to the total spectrum and, hence, it is not 
easy to extract the contributions of specific configurations any more. 
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It is worth pointing out here that the calculated results of the total 
spectra in Fig. 10 are similar for all cases shown here insensitively to 
the photon energy. Thus, we need to keep in mind that the contribu- 
tion of each subcomponent has certain energy dependence even if the 
total spectrum has almost no energy dependence. 

Next let us show the (7,p) spectra calculated with the attractive 
potential V^ 1 and the repulsive potential vif 1 defined in Eqs. (9) (10). 
We evaluate the spectra with these potentials at the incident photon 
energy = 2.7 GeV. As seen in Fig. 11(a), in the case of the attrac- 
tive potential vS°^ , the bound state structure is clearly seen thanks to 
the narrow width of the bound state in p state and dominant contri- 
bution of the (0p 3 / 2 )p 1 <£>Pu> configuration. This is consistent with the 
(7,p) spectrum shown in Ref. [12], in which the energy dependence on 
the oj optical potential is taken into account. Comparing these two 
spectra, we see that the energy dependence of the potential has rela- 
tively small contribution to fix the global shape of the spectra and the 
energy dependence affects on the reaction rate as a weak repulsion. 
In Fig. 11(b), we show the (7,p) spectrum with the repulsive poten- 
tial. Apparently the shape of the spectrum is quite different from the 
previous two attractive potentials. The repulsion of the omega meson 
inside the nucleus makes the spectrum push to higher energies and 
almost nothing is seen below the omega threshold. It would be in- 
teresting if one could measure the omega-nucleus system in the (7,p) 
reaction around the threshold of the omega production. 

Here we have obtained three types of plots for the formation spec- 
tra of the omega-nucleus system. In the first case, the optical potential 

(a) 

Vu} has enough attraction to form bound states of the u> meson in nu- 
cleus, but the widths of the bound states are enlarged due to the strong 
absorption. Thus the whole spectrum of the (7,p) reaction dose not 
show distinct peak structure caused by the bound states even at the 
ideal recoilless kinematics, in which the {^Pz/2)p l ® Pu> configuration 
dominants the spectrum. In the second case, since the optical poten- 
tial vS°^ has more attraction and less absorption, a deeper bound state 
of omega in p wave is formed with a narrower width. This bound state 
is clearly seen in the spectrum at the recoilless condition. In the third 

(c) 

case, the optical potential is repulsive. The spectrum obtained 
with this potential has definitely different from the above two cases. 
The missing mass spectroscopy of the (7,p) reaction in the omega pro- 
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duction energies is a good practical tool for the investigation of the 
omega-nucleus interaction, at least it will be clearly seen if the optical 
potential is attractive or repulsive. 

4.3 cr-nucleus system 

We propose the (7,p) reaction as a new method to create the sigma 
meson inside a nucleus, as discussed in Ref. [18] with the (d, 3 He) 
reaction. In this subsection, we discuss the sigma meson production in 
the (j,p) reaction on a heavy nucleus target 208 Pb. We show here the 
total spectral function formulated in Sec. 3.2 instead of the differential 
cross section because of lack of our knowledge for the elementary cross 
section of the sigma meson photoproduction. Since the total spectral 
function involves the effective number which counts the reaction rate 
of the picked-up proton induced by photon, we can investigate an 
impact of the narrow bound states of the sigma meson around the 
two pion threshold, showing the global shape of the spectral function 
without a definite number of the cross section. The incident photon 
energy in the (7,p) reaction is chosen to be 400 MeV so as to satisfy the 
recoilless condition for two-pion production, since we are interested in 
the deeply bound sigma meson. 

Fist of all, we show the total spectral functions p^ot ( w ) f° r the 
(d, 3 He) reactions with the 208 Pb target for comparison. The incident 
deuteron energy is chosen to be 1.5 GeV in order to satisfy the re- 
coilless condition around the emitted particle energy Tf = T{ — 2m n . 

Figure 12(a) shows the plots of p to [ ' {uj) for the C = 0.2 and 0.4 
cases. In the (d, 3 He) reaction, one finds a less prominent peak around 
T{ — Tf = Im-x for C = 0.4 as reported in Ref. [18]. It would be hard, 
however, to distinguish such a small peak in experiment. One sees also 
that the whole shape of the spectral function for C = 0.4 resembles 
that for C = 0.2. This is because the large distortion effect of the 
hadronic probe suppresses strongly production of the deeply bound 
states, which have the significant enhancement around u> ~ 2m n , and, 
hence, the total spectral function is dominated by the shallow bound 
states with large widths like the C = 0.2 case. 

Next, in Fig. 12(b), we show the total spectral functions for the 
(7,p) reaction on the 208 Pb target with C = 0.2, 0.3 and 0.4 as a 
function of the missing mass Tj — Tf. As shown in the figure, the 
spectral function of the C = 0.4 case is significantly enhanced around 
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2m n in a reflection of production of the deeply bound states thanks 
to the transparency of the (7,p) reaction. The incident photon can go 
into the center of the nucleus, knocking out a proton lying in a deep 
state, and, as a consequence, the deeply bound sigma states with 
narrow widths are formed in the nucleus. This is the good advantage 
of the photon induced reaction. The figure shows also a peak structure 
of the spectral function with C = 0.3 above the two-pion threshold. 
This is also the consequence of production of the deeply bound sigma 
states. 

In order to compare the energy dependence of the total spectral 
functions for the different strength of the chiral restoration, which is 
parameterized as C, we show in Fig. 12(c) spectral functions normal- 
ized so as to have the almost same peak height as the C = 0.4 case. 
The normalization factors are 10 for C = 0.3, 40 for C = 0.2 and 120 
for C = 0.0. For comparison, we show also a spectral function with 
C = 0.0, in which the partial restoration of chiral symmetry does not 
take place. The spectral function with C = 0.0 is obtained in Ref. [27], 
having the peak around Tj — Tf ~ 550 MeV. Note that the spectral 
function with C = 0.0 is calculated without the effective number and, 
therefore, the reaction mechanism is not involved in the spectral func- 
tion unlike the cases of C ^ 0.0. Nevertheless the comparison gives 
us a schematic view of the energy dependence. 

We can see the clear shift of the peak energy in this figure and can 
expect that the shift could be a good signal of the chiral restoration 
in the nuclear medium. 

4.4 Background Consideration 

So far we discuss the production rate of the mesons in nucleus in 
the (7,p) reaction. In the missing mass spectroscopy of the emitting 
proton, the production of the meson in the nucleus is expected to be 
seen on top of a certain structureless background. In this section, we 
roughly estimate the background of the (7,p) reactions. 

First of all, in the ?7-nuclear system, the typical signal for the cre- 
ation of the eta meson in nucleus is estimated to be 100 nb/(sr MeV) 
for the differential cross section of the 12 C(7,p) reaction at the recoil- 
less energy = 950 MeV as discussed in Sec. 4.1. On the other 
hand, experimental data [51] showed that the differential cross section 
of the inclusive 12 C(7,p) reaction is d 2 a/d£ldp ^ 0.9 [/xb/(sr MeV/c)] 
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at the photon energy E y = 580 MeV and the emitted photon angle 
dp = 23° in the laboratory frame. From the data, we estimate the 
physical background due to the inclusive proton emission is 



E d 2 



a 



1.5 



fib 



(27) 



dQdE p dVtdp 



sr MeV 



in the proton energy distribution. Therefore, the signal over noise 
ratio is expected to be around 1/15. However, the kinematical con- 
dition of the (7,p) reaction in Ref. [51] is different from that of the 
7/-nucleus system formation reaction in Section 4.1. Hence, we need 
to use theoretical models and/or more appropriate experimental data 
to improve the background estimation. 

As for the u> mesic nuclei formation by the (7,p) reaction, the sig- 
nal is estimated to be 10 nb/(sr MeV) for the 12 C target case in 
Sec. 4. 2, while the background is expected to be smaller than 100 
nb/sr/MeV [24]. Thus, the signal over noise ratio is to be around 
1/10 from the calculated results. 

For the a- nucleus systems, we do not have any background es- 
timations so far, since we cannot evaluate definite values for the a 
meson production rates due to lack of the elementary cross section. 
Conceivable physical backgrounds would be resonant enhancements 
of the cross section due to nucleon excitations in nucleus. Especially 
the kinematics considered here of the two pion production could con- 
flict with a A excitation in nucleus. A possible way to avoid such a 
physical background is to observe more particles in the final state, for 
instance, with taking coincidence of the emitted proton and two pion 
with isospin 0. Another possibility might be to observe the (7,p) re- 
actions with lighter nuclear targets and to see the mass dependence of 
the peak around the two-pion threshold, since the deeply bound sigma 
states will not be seen in a lighter nucleus. The detailed discussions of 
the physical backgrounds are beyond the scope of the present paper. 



In this paper we have made a theoretical evaluation of the formation 
rates of the rf and to mesons in nuclei induced by the (7,p) reactions 
in ideal recoilless kinematics. We have shown the expected spectra in 
order to investigate the meson-nucleus interactions. We have found 



5 Conclusion 
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that the (7,p) reactions are good practical tool to investigate the prop- 
erties of the mesons created deeply inside the nucleus due to the small 
distortion effects. This good advantage provides the distinct differ- 
ence in the formation spectra of the 77-nucleus system obtained by the 
two chiral models which are based on the different physical pictures of 
the iV(1535) resonance. For the w-nucleus system, we have compared 
three types of the 00 optical potentials in the (7,p) spectra, showing 
the definitely different shapes of the spectra. 

We have also investigate the (7,p) spectra at the recoilless con- 
dition for the two pion production in isoscalar channel in order to 
study an impact of the creation of the deeply bound states of the 
sigma meson associated with the partial restoration of chiral symme- 
try in heavy nuclei. We have found that a prominent enhancement 
around the two pion threshold in the missing mass spectra in case of 
a sufficient strength of the partial restoration in medium, owing to 
the transparency of the (7,p) reaction to create deeply bound states. 
Thus we expect that the (7,p) reaction is a good tool to create the 
sigma meson in nucleus. 

The study of the bound states is one of the most promising method 
to investigate the meson properties at finite density. Nevertheless the 
large natural widths of the meson bound states in nucleus disable 
to distinguish contributions from each bound state. In such 
global conformation of the missing mass spectra is necessary to ex- 
tract valuable information of the meson nucleus interaction. It is also 
beneficial to compare the spectra of the (7,p) and (d, 3 He) reactions 
at corresponding recoilless conditions, since each configuration differ- 
ently contributes to the total spectra due to the different distortion 
effects. We expect that the present results stimulate the experimental 
activities and help the developments of this research field. 
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Figure 8: Calculated spectra of 12 C(d, 3 He) n B®?7 reactions at Td = 3.5 GeV 
in the left-side and 12 C(7,p) n B®?7 reactions at E 1 = 950 MeV in the right- 
side are shown as functions of the excitation energy E ex defined in the text. 
Eq is the rj production threshold energy. The 7/-nucleus interactions are 
evaluated by (a) the t-p approximation, (b) the chiral doublet model with 
C = 0.2 and (c) the chiral unitary model [15]. The total spectra are shown by 
the thick solid lines, and the contributions from the dominant configurations 
(0si/2) p 1 <8> s v and (0p 3 / 2 )p 1 ®p n are shown by dashed lines and dash-dotted 
lines, respectively, where the proton-hole states are indicated as {nij)' 1 and 
the rj states as t n . 
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Figure 9: Calculated spectra for the formation of the r^-mesic nuclei in (a) 
the (d, 3 He) reaction at T d = 3.5 GeV and (b) the (7,p) reaction at E y = 950 
MeV are shown as functions of the excitation energy E ex defined in the 
text. Eq is the rj production threshold energy. The 7/-nucleus interaction 
are evaluated by the chiral doublet model with C = 0.0 and C = 0.2, and 
the expected spectra are shown by thin lines and thick lines, respectively. 
The total spectra are shown by the solid lines and dominant subcomponents 
are shown by dashed lines as indicated in the figures. Here, the proton-hole 
states are indicated as {nij)' 1 and the rj states as £ v . 
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Figure 10: Calculated spectra of 12 C(7,p) reaction for the formation of the 
o n B system with the attractive potential at (a) E 1 = 1500 MeV, (b) 
£ 7 = 2000 MeV, (c) £ 7 = 2700 MeV and (d) £ 7 = 3000 MeV as functions 
of the excitation energy E ex defined in the text. E is the uj production 
threshold energy. The thick solid lines represent the total spectra and the 
other lines represent the contributions from the dominant subcomponents as 
indicated in the figures. 
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Figure 11: Formation spectra of the u mesic nucleus in the (7,p) reaction 
with the 12 C target calculated with the attractive potential (a) and 
the repulsive potential (b) at the photon energy E 1 = 2.7 GeV. E is 
the ijj production threshold energy. The thick solid lines represent the total 
spectra and the other lines represent the contributions from the dominant 
subcomponents as indicated in the figures. 




Figure 12: Total spectral functions ptot defined in the text for (a) (d, 3 He) re- 
actions and (b) (7,p) reactions on 208 Pb target. The incident particle energies 
Tj are (a) 1.5 GeV and (b) 400 MeV, respectively. The solid lines indicate 
the spectral functions with the parameter C = 0.4. The dashed and the 
dash-dotted lines represent the spectral functions with C = 0.3 and C = 0.2, 
respectively, (c) Normalized spectral functions of the (7,p) reactions. The 
normalization factors are 10, 40, 120 for the C = 0.3, C = 0.2 and C = 0.0 
cases, respectively. The spectral function with C = 0.0 reported Ref. [27] is 
also shown for comparison. 
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